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ABSTRACT 

Applying our recently-developed generalized version of the high-latitude emission theory to the 
observations of X-ray flares in gamma-ray bursts (GRBs), we present here clear observational evidence 
that the X-ray flare emission region is undergoing rapid bulk acceleration as the photons are emitted. 
We show that both the observed X-ray flare light curves and the photon index evolution curves can 
be simultaneously reproduced within a simple physical model invoking synchrotron radiation in an 
accelerating emission region far from the GRB central engine. Such an acceleration process demands 
an additional energy dissipation source other than kinetic energy, which points towards a significant 
Poynting-flux in the emission region of X-ray flares. As the X-ray flares are believed to share a similar 
physical mechanism as the GRB prompt emission, our finding here hints that the GRB prompt 
emission jets may also carry a significant Poynting-flux in their emitting region. 

Subject headings: gamma-ray burst: general — radiation mechanisms: non-thermal — relativistic 
processes 


1. INTRODUCTION 

Gamma-ray bursts (GRBs), the most luminous explo¬ 
sions in the universe, invoke relativistic jets beaming to¬ 
wards Earth with the highest velocities for b ulk motion 
in the universe (e.g., Kumar &; ZhandI2015H . Some of 
them are followed by softer, less energetic, X-ray flares, 
which also move with relativistic velocities towards Earth 
(jBurrows et al.ll2005t ). The X-ray flares are detected by 
the Swift satellite at ~ (10 2 — 10 5 ) seconds after the 
GRB trigger, show a rapid rise and a steep fall, out¬ 
shine the underlying afterglow by up to ~ 3 orders of 
magnitude in energy flux, and display a distinctive pat¬ 
tern of spectral evolution with a strong spectral soften¬ 
ing during the decay phase. All these features are in¬ 
consistent with the standard afterglow emission, emitted 
from an external blas t wave that propagates through a 


surrounding medium (Meszaros & Rees 19971: Sari et al. 

1993). Indeed, observations (e.g., 

Liang et al.l 2006); 

Chincarini et al. 20071:IMarerutti et al.l 

201® and theoret- 

ical modelings (e.g., 

Zhang et ajj 2006; Lazzati & Perna 


l2007t iMaxham fe Z hang 2009) rather suggest that the 


X-ray flares share a similar physical mechanism as the 
GRB prompt emission itself, which results from internal 
energy dissipation within the jet. 

For a spherical, relativistic jet, the decay of light curve 
cannot be steeper than a certain level defined by the 
so-called “high-latitude emission effect”, or the “curva¬ 
ture effect”. This is to say, even if the emission stops 
abruptly, photons from higher latitudes with respect to 
the observer’s line of sight would arrive the detector at 
progressively later epochs but with a progressively lower 
Doppler factor value, resulting in a steep decay in flux. 
If the emission region keeps a co nstant L orentz factor 
T, there exists a simple relation (iKurnar fe Panaitescul 

Earn 

a = 2 + /3 (1) 

between the temporal index a and the spectral index 
/?, in the convention that the observed spectral flux is 


expressed as oc t “ b “ v ~^ s , where t 0 bs is the observer 

time and u 0 b s is the observed frequency. Such a curvature 
effect has been invoked to int erpret th e fas t decline of the 
X-ray flare light curves dLiang et alJl2006f h 

In this Letter, we confront the steep decay phase of 
X-ray flares with our g eneralized version of the high- 
latitude emission theory (lUhm fe Zh ang 2015) and show 
that the emitting region of X-ray flares undergoes rapid 
bulk acceleration. Also, for the first time, we present a 
physical modeling to both the observed flare light curves 
and the photon index evolution curves, simultaneously. 

2. HIGH-LATITUDE CURVATURE EFFECT AND THE 
STEEP DECAY PHASE OF X-RAY FLARES 

In practice, testing the curvature effect theory with 
data is obscured by the so-called “T 0 -effect” flZhang et al.l 
120061 1. Since GRB light curves are plotted in logarithmic 
scale for both the observer time and the flux, the appar¬ 
ent decay slope dobs sensitively depends on the reference 
time (t 0 bs = 0) to plot the light curves. Usually t D b s = 0 
is defined at the GRB trigger time. For an X-ray flare, 
however, the emission episode likely starts at a later time 
(say, t 0 bs = To) than prompt emission, so that a new zero 
time (t 0 bs = 0 at To) should be chosen in order to cor¬ 
rectly study the radiation physics of the flare. Keeping 
the GRB trigger time as time zero point would cause an 
un-physically rapid decay that is steeper than the the¬ 
oretical prediction. Indeed, the observed photon index 
T in the XRT band during the decay phase of an X-ray 
flare is typically in the range of ~ (1 — 3), so that the 
corresponding spectral index, /3 = T — 1, is in the range 
~ (0 — 2). The expected temporal index d from Equation 
m should then be in the range ~ (2 — 4), yet the ob¬ 
served decay indices d 0 bs are usually significantly larger 
than these values (up to more than 10). This has been at¬ 
tributed to the T 0 effect, suggesting that the X-ray flares 
have l ater emissio n epis odes than the GRB prompt emis¬ 
sion (|Liang et al.ll2006h . 

Correctly accounting for the To effect is not straight¬ 
forward, because the true beginning of the flares cannot 
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be identified directly from the data. Depending on the 
detector sensitivity, the actual beginning may have been 
missed since it is too faint to be detected. More impor¬ 
tantly, the initial portion of the flares is usually buried 
under the background afterglow emission. Thus, one may 
formally write that the beginning of an X-ray flare is at 

T 0 = T 0 obs - (AT) missed , (2) 

where T^ hs is the “observed” starting time of the flare 
and (AT) mlssed is the time duration of the missed por¬ 
tion. If one defines a new zero time fobs = 0 at To, the 
new light curve of the flare would exhibit a shallower 
decay than the original observed one. The larger the 
T 0 value one corrects for, the shallower the decay slope 
one would get. The true To value cannot be constrained 
observationally due to the unknown quantity (AT) mlssed 
but may be determined through theoretical modeling. 

We analyze three X-ray flares observed in GRB 
140108A, GRB 110820A, and GRB 090621A (see Fig¬ 
ure [[]). The data are taken directly from the UK 
Swi ft Science Data Ce ntre at the University of Leices¬ 
ter (lEvans et al.ll2007T) . and the 10 keV F light curves 
(top panels) and XRT-band (0.3-10 keV) photon index 
curves (bottom panels) are presented. By identifying a 
flare (blue data points) from an original light curve of 
each example, we define Tp bs as the first data point of 
the flare and mark it with the dotted vertical line in the 
top panel. We first take T 0 equal to Tp bs , which is the 
upper limit value of To. By systematically subtracting 
To from f 0 b s of each data point of the flare, we effectively 
re-define the zero time point and “translate backwards” 
the data points to construct a new light curve of the X- 
ray flare (red data points). The observed photon index 
(T) curve (blue data points) during the identified flare 
region and its corresponding backward translation (red 
data points) are also presented. 

By fitting the photon index curve during the decay 
phase of the flare (green curve in the bottom panel), we 
reconstruct the predicted decay light curve (green curve 
in the top panel) based on Equation (fTj). which has the 
steepest decay slope allowed by the curvature effect as¬ 
suming that the jet is moving with a constant speed. One 
can see that the observed decay light curve violates this 
limit significantly in all three examples. Strictly speak¬ 
ing, a should be derived by fi at 10 keV (since the light 
curves are at 10 keV). The measured photon indices, and 
hence f3, on the other hand, are defined for the XRT 
band. A possible correction may be made through de¬ 
tailed modeling. Using the modeling results for the GRB 
140108A X-ray flare (Figure [2] for details) as an exam¬ 
ple, we calculate the 10 keV time-dependent photon in¬ 
dex (black curve in the bottom panel of GRB 140108A). 
Based on Equation 0, the predicted decay light curve is 
constructed again (black curve in the top panel of GRB 
140108A). Even though it is slightly steeper than the 
green curve, it is still significantly shallower than the ob¬ 
served light curve. Performing the same exercise to the 
other two X-ray flares leads to the same conclusion. 

Notice that we have adopted the most conservative 
value for To (i.e., its upper limit, Tg bs ) in performing the 
test whether the observed decay phase of X-ray flares is 


consistent with Equation 00 Correcting for a more re¬ 
alistic value (smaller than Tg bs ) would lead to an even 
steeper decay curve after backward translation, worsen¬ 
ing the conflict with Equation 0 discussed above. 

In a recen t work on high-latitude emission theory 
(lUhm fc Zhandr2015ll . we find that Equation 0 strictly 
holds only for a constant value of the Lorentz factor T. 
In the case of an accelerating shell, the temporal decay 
index a becomes significantly larger than 2 + (3. The 
trend is opposite for the case of a decelerating shell. Ap¬ 
plying this new generalized version of the high-latitude 
emission theory to the apparent tension between the ob¬ 
servational data and Equation 0, we immediately con¬ 
clude that the emission regions of the three X-ray flares 
are undergoing significant bulk acceleration as the X-ray 
photons are emitted. 


3. A SIMPLE PHYSICAL MODEL FOR X-RAY FLARES 

In order to verify the statement above, we perform 
detailed numerical modeling to these three X-ray flares 
shown in Figure |T] We aim at reproducing the entire 
temporal behavior, i.e., not only the decaying phase but 
also the rising phase, for both the flare light curve and 
the photon index curve. The redshifts of the three GRBs 
are all unknown, and we assume a typical value z = 1. 
The redshift plays o nly a global role in s haping the ob¬ 
served spectral flux (Uhm fc Zhangll2015l) . and therefore 
this assumption would not degrade the outcome of our 
modeling. The standard flat ACDM universe with the 
parameters Hq = 71 km s” 1 Mpc -1 , D m = 0.27, and 
Da = 0.73 is adopted in our calculations. 

We adopt a simple physical picture: a single relativistic 
spherical shell expands radially with a profile of the bulk 
Lorentz factor T(r) as a function of radius r. The shell 
starts to emit photons at radius r on (and at the lab-frame 
time ton), and finishes emitting at radius r 0 g. During the 
emission phase, photons are continuously emitted from 
all locations in the shell, with an isotropic angular distri¬ 
bution of the emitted power in its co-moving fluid frame. 
Considering synchrotron radiation as the radiation mech¬ 
anism (as sug gested by recent theoretical modeling an d 
data analysis (lUhm fc Zhand 120141 : iZhang et all 1201 bil l, 
we delineate the shape of the ph oton spectr um in the 
fluid frame in the form of ijtJhm X Zhanell20Tfljl 

H(x) with x = v' jv ' h , (3) 


with the comoving-frame characteristic synchrotron fre¬ 
quency 


y ch 


_3_ q e B 2 

16 m e c ^ ch ’ 


( 4 ) 


where m e and q e are the mass and charge of the electron, 
respectively, and c is the speed of light. The magnetic 
field strength B and the characteristic Lorentz factor 


1 In case where the photons are emitted from a region with a 
finite width, Tq° s gives the upper limit value of To only for the first 
shell in the emitting region while the test needs to be performed 
for the last shell. However, if an X-ray flare is produced from such 
a shell with a finite width, an expected feature of turning-off of the 
shell emission would be a rather flat shape of the peak area, due to 
the contributions from all shells in the region. All three examples 
in Figure [l] show a relatively sharp peak, without indicating such 
a flat feature at the peak area. Therefore, the analysis presented 
here remains valid. 
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7ch of the electrons are measured in the fluid frame, co¬ 
moving with the shell. The observed spectral flux, Fj,° bs , 
is calculated based on Doppler transformations from the 
comoving frame to the observer fr ame, with the curva - 
ture effect fully taken into account (lUhm fc Zhandl20T5lf . 
The number of radiating electrons in the shell is assumed 
to increase at an injection rate i?i n j (measured in the fluid 
frame) from an initial value N = 0. 

If one resets the reference time at To, the observer- 
frame time t 0 bs of the flare can be physically connected 
to the lab-frame time t through 


bibs — [bin T c(t t on ) T COS (1 f Z ) 

( t_ c M ) _ (* on ~1r)] (1+2)i 


(5) 


taking into account the point that the initial photons of 
the flare are not emitted at r = 0, but rather at r on at 
ton- Here, fi = cos 0, and 6 is the polar angle measuring 
the latitude of emission location with respect to the ob¬ 
server’s line of sight. Notice that the observer time fobs 
is, in fact, independent of t on since the time t is calcu¬ 
lated as t = ton + f r dr/(cl 3), where /3 = (1 — 1/T 2 ) 1 / 2 . 

However, in practice, the true To value cannot be con¬ 
strained directly from the data. Hence, we may instead 
reset the reference time at Tg bs , i.e., at the observed 
starting time of the flare. Then, according to Equation 
©, the observer time of the flare should read 


^obs — 





(1 + z)~ (AT) missed . 


(6) 

This observer time is supposed to describe the red data 
points in Figure [1] since those red points were obtained 
by reseting the reference time at Tg bs . Thus, we use 
Equation © and model the red data points theoretically, 
which would enable us to constrain the unknown missed 
portion (AT) missed . 

In our modeling, we introduce power-law dependences 
of the following parameters: 


T(r’) =T 0 (~) > (7) 

BM=5oW , (8) 

7ch(r)=7ch(^) ■ (9) 


The index s delineates the degree of acceleration. The 
index 6, with a typical value of 1, is naturally expected 
due to H-field flux conservation in an expanding shell, 
which is crucial to interpret GRB spectrum in the fast - 
cooling regime (jllhm fc Zhandl20l4 IZhang et al.ll201fiH . 
The index g describes how the characteristic electron 
Lorentz factor evolves with radius r. The shape of 
the photon spectrum H(x) may be taken as a cutoff 
powe r law, or the so-called “Band” function (jBand et all 
Il993fl. which usual ly well describes the spectra of GRBs 
( Preece et al.ll2000l f and can be accounte d for within the 
fast-cooling synchrotron radiation model (Uhm & Zhang 
[Ml IZhang eiTaD[20ill . 

We begin our calculations at the radius r on . By setting 
r 0 = r on and taking 6=1 and Bq = 300 G, we search for 


s, g , r 0 ff, i?inj, H(a ;), and (AT) mlssed to reproduce the 
observation of each of the three X-ray flares. The results 
are shown in Figure [2j In all three examples, one can see 
a good agreement between our modeling and the data for 
both the lOkeV-Fb light curves (top panels) and the pho¬ 
ton index evolution curves (bottom panels). The model- 
predicted light curves at 1 keV (green) and 3 keV (cyan) 
are also presented in the top panels. A time-dependent 
model-prediction curve for the XRT-band photon index 
is calculated as a power-law index between 1 keV and 10 
keV light curves (black curve in the bottom panels). We 
note that without introducing the acceleration index s, 
the light curves and the photon index curves could not 
be reproduced. This suggests that our conclusion of an 
accelerating emission region for the X-ray flares in GRB 
140108A, GRB 110820A, and GRB 090621A is robust. 

Following parameters are adopted in common in all 
three examples: r on = ro = 10 14 cm, B 0 = 300 G, and 
6 = 1. The spectral function is taken as a power law 
with an exponential cutoff, i.e., H(x ) = x a+1 e~ x with 
a = —0.7. A Band-function (lBand et all I1993T ) input 
spectrum with a steep high-energy index gives a simi¬ 
lar overall match to the data. Other model parameters 
for each example are as follows. For (GRB 140108A, 
GRB 110820A, GRB 090621A): r 0 = (7.5,6.0,7.0), 
s = (1.15,0.95,0.8), 7 ° h = (0.82,2.2,2.2) x 10 4 , g = 
(0.48,0.0,0.05), r-off = (2.0,2.0,1.5) x 10 15 cm, and a 
constant injection rate, T2 in j = (6.0,2.5,5.5) x 10 47 s _1 . 
A non-zero (AT) mlssed is indeed required to model the 
data, and for each case, (AT) mlssed = (18,12,36) s is 
adopted. 

The rising phase of the flare in GRB 140108A exhibits 
a strong spectral hardening; the photon index decreases 
from F ~ 2.4 to f ~ 1. This implies that the peak energy 
of the observed photon spectrum rapidly increases during 
this phase. According to our modeling, the spherical 
shell is actively emitting during the rising phase, and the 
temporal behavior of the observed peak energy roughly 
follows E° hs ~ TB^ 2 ch ~ r s - fc+2g . For the flare in GRB 
140108A, the result of our modeling gives E° hs ~ r 111 , 
confirming the rapid evolution of the peak energy during 
the rising phase. On the other hand, the strong spectral 
softening observed during the decay phase of all three 
flares is due to the high-latitude curvatu re effect in the 
acceleration regime (lUhm fo ZhandlMbil . The emission 
from higher latitudes has a progressively smaller Doppler 
boosting, and therefore the observed spectrum is placed 
with a progressively smaller peak energy, resulting in a 
spectral softening. 

4. CONCLUSIONS AND DISCUSSION 

In this Letter, by applying our newly-developed gen¬ 
eralized version of the high-latitude emission theory 
{Uh m fe Zhanell2015( ) to the observations of three exam¬ 
ple X-ray flares in GRB 140108A, GRB 110820A, and 
GRB 090621A, we have presented clear observational 
evidence that the emitting region of these X-ray flares 
should undergo rapid bulk acceleration as the flare pho¬ 
tons are emitted. Furthermore, we have shown that the 
entire observed temporal behavior, i.e., not only the de¬ 
caying phase but also the rising phase, for both the flare 
light curves and the photon index evolution curves can 
be simultaneously reproduced within a simple physical 
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model invoking synchrotron radiation in an accelerating 
emission region far (~ 10 15 cm) from the GRB central 
engine. This is the first time that such a comprehensive 
theoretical modeling is done for GRB X-ray flares. 

The identification of an acceleration process in a rel¬ 
ativistic jet has profound implications. One important 
opening que stion in the field o f GRBs is the composition 
of the jets (jKumar fc Zhanj I2015H . Within the stan¬ 
dard matter-dominated “fireball” m ode l, t he jet un der¬ 
goes rapid acceleration early on (IMeszaros et al.l 1 1 9931 
iPiran et ~ail 1 9931 iKobavashi et al.lll999H below the coast¬ 

ing radius r c < Ti?o = 3 x 10 11 cm(T/300)(i?o/10 9 cm) 
(where Rq is the radius of the jet base), and becomes 
kinetic energy dominated afterwards. The non-thermal 
emission is believed to be emitted from internal shocks 
wher e the kinetic energy is dissipated (|Rees fc Meszarosl 
Il994f) . Within this scenario, the emission is released 
at the expense of the kinetic energy, so that after the 
prompt emission phase, the average Lorentz factor of the 
fireball is expected to be reduced. An alternative sce¬ 
nario invokes a Poynting-flux-dominated outflow, with 
the magnetization parameter a (the ratio between Poynt- 
ing flux and matter flux) greater than unity in the emis¬ 
sion region. It has been known that such kind of jet may 
undergo slow acceleration eve n without magnetic dis¬ 
sipation and photo n radiation (jKomissarov et all 120091 
iGranot et al.l 1201 1 ), as the Poynting flux energy is grad¬ 
ually converted to kinetic energy due to the magnetic 
pressure gradient within the jet. If the Poynting flux un¬ 
dergoes an abrupt dissipation, probably due to internal 
collision-ind uced magnetic reconnection and turb ulence 
(ICMART) dZhang fc Yanll2f)ITt iDeng et al.ll2015H . part 
of the dissipated Poynting flux energy would be given 
to the jet for bulk acceleration as the other portion of 
the energy is converted to particle energy and released 
as photons. Such a model also invokes a relatively l arge 
emission radius from the central engine (jZhang fc Yanl 
l 201 lh . which is required to interpret the long-duration 
X-ray fl are d ecay t ail within the curvature effect model 
dZhang et al.ll2006h . As a result, our finding provides a 
“smoking-gun” signature of a > 1 in the emission region 


of X-ray flares. This aligns with other argume nts that X - 
ray flares are Poynting-flux-dominated dFan et al.ll2005f) . 


Observational and theoretical arguments suggest that 
GRB prompt emission shares a si milar p hysical ori- 


gin with the X-rav flare emission ( 

Liang et al. 2006t 

Chincarini et al. 2007; Margutti et alj 

2010; Zhang et alj 

2006| Lazzati & Perna 20071; Maxham & Zhang 2009)). 


Therefore, our finding of bulk acceleration in X-ray 
flares hints that the emission region of GRBs may also 
have er > 1 , which is consistent with other observa¬ 
tional evidence, including the lack of or the weak t her¬ 
mal emission component in most GRBs (Zh ang fc Pe'erl 
120091: fGao fc Zhanel [2015ft . the polarized 7 -ray emission 
( Yonetoku et all I2011T) . and the polarized early optical 
afterglow emission from t he reverse shoclfl (' Steele et all 
l2009t Mundell et al.ll201~3h . Indeed, based on an indepen¬ 
dent argument regarding the observed spectral lags, we 
recently showed that the GRB prom pt-emissio n region is 
also undergoing bulk acceleration (lUhm fc ZhanetramTi) . 

A decay index a steeper than 2+/3 can also be achieved 
by inv oking anisotropic emissio n in the jet co-moving 
frame dBeloborodov et aD 12011 ). However, both spec¬ 
tral and temporal properties of X-ray flares, i.e., strong 
spectral evolution observed during both the rising and 
decaying phases of flares as well as steep decay light 
curves, cannot be easily interpreted within such a sce¬ 
nario. Bulk acceleration is demanded in order to simulta¬ 
neously interpret the light curve an d spectr al evolution of 
the flares. An extended analysis (j.Tia et al.ll2015h shows 
that the majority of X-ray flares are in the acceleration 
regime, suggesting a ubiquitous Poynting-flux-dominated 
composition among X-ray flares. 
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Fig. 1.— Three example X-ray flares in GRB 140108A, GRB 110820A, and GRB 090621A, that show strong observational evidence of 
significant bulk acceleration in the X-ray flare emission region. In the top panels, the black and blue data show the 10 keV spectral flux 
light curve as observed by XRT, with the blue color being used to identify the X-ray flare part. The blue data in the bottom panels show 
the photon index of the XRT band during the identified flare region. The dotted vertical line in the top panels represents the beginning of 
the observed data for each X-ray flare, i.e., Tg hs . We translate blue data points backwards by the amount Tg hs and obtain red data points 
in top and bottom panels, respectively. By fitting the photon index evolution during the decay phase of flares (green curve in the bottom 
panels), we construct the predicted decay light curve (arbitrary normalization) based on Equation {TJ and show it by the green curve in 
the top panels. A more stringent approach may be possible through detailed modeling, as indicated by the black curves in GRB 140108A 
X-ray flare. See Se ction [2] for details. All the original data are taken from the UK Swift Science Data Centre at the University of Leicester 
(|Evans et alJl2007n . 



Fig. 2.— Modeling the three X-ray flares in GRB 140108A, GRB 110820A, and GRB 090621A, using the curvature effect with an 
accelerating emission region. Top panels present the observed flare light curve at 10 keV (red data points) and our model light curve at 10 
keV (black curve). The model-predicted light curves at 1 keV (green) and 3 keV (cyan) are also presented. Bottom panels present a model 
prediction to the XRT-band photon index (black curve) and the observational data (red data points). See Section [3] for details of modeling. 
The blue data points indicate the original observed light curve (top panel) and photon index curve (bottom panel) of each X-ray flare. 








































